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Abstract 

A number of two-dimensional(2D) critical phenomena can be described 
in terms of the 2D sine-Gordon model. With the bosonization, several 
ID quantum systems are also transformed to the same model. How- 
ever, the transition of the 2D sine-Gordon model, Berezinskii-Kosterlitz- 
Thouless(BKT) transition, is essentially different from the second-order 
transition. The divergence of the correlation length is more rapid than 
any power-law, and there are logarithmic corrections. These pathological 
features make difficult to determine the BKT transition point and critical 
indices from finite-size calculations. In this paper, we calculate the sev- 
eral correlation functions of this model using a real-space renormalization 
technique. It is found that the several correlation functions, or eigenvalues 
of the corresponding transfer matrix for a finite system, become degen- 
erate on the BKT line including logarithmic corrections. By the use of 
this degeneracy, which reflects the hidden SU(2) symmetry on the BKT 
line,it is possible to determine the BKT critical line with high precision 
from small size data, and to identify the universality class. In addition, 
a new universal relation is found. This reveals the relation between the 
Abelian and the non-Abelian bosonization. 

PACS numbers: 75.40.Cx, 05.70.Fh, ll.lO.Hi, 75.10.Jm 
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1 Introduction 



Sine-Gordon model in two-dimension(2D), which is a natural extension of the 
Gaussian or free boson model, plays important role in 2D classical and ID 
quantum systems, such as 2D XY model, 2D Helium film, ID quantum spin 
models and fermion models. 

The peculiarity of the phase transition of the 2D XY model(Helium filni) are 
first noticed with the spin wave approximation or the free boson model||l|, ^. 
In these theory, there is no continuous symmetry breaking as expected by the 
Mermin- Wagner theorem||^, but at all temperatures the correlation length ^ 
is infinite and the correlation functions decay as power-law whose exponents 
vary continuously with temperature. This picture is qualitatively correct at low 
temperatures, but clearly wrong in high temperatures where one expects a finite 
^ and the associated exponential decay of correlations. 

BerezinskiiQ and Kosterlitz and Thouless[^ pointed out the importance of 
vortex excitations, which modify essentially the spin wave theories. The vor- 
tex structure reflects the periodic nature of the spin variable (f) = (f> + 2?i7r. 
The vortices carry integer vorticity and interact among themselves via the log- 
arithmic 2D Coulomb interaction. At low temperatures, all of the particles are 
bound into neutral "quasi-molecules" with zero vorticity, so they only change 
the coupling of the spin wave model. At higher temperatures, the binding of 
the quasi-molecules decreases and it causes a phase transition. 

Kosterlitz]^ subsequently performed a renormalization group calculation, 
following the method by Anderson et al.|0, ||]. In fact, the structure of the 
renormalization equations is the same for both cases. He found that close to the 
Berezinskii-Kosterlitz-Thouless(BKT) transition point, the correlation length 
diverges as 

C cx exp(5Vi), (1) 

faster than any power of t. Also on the BKT critical point there appear loga- 
rithmic corrections in various quantities, such as correlation functions and sus- 
ceptibilities. These features are entirely different from the usual second-order 
transition. 

The 2D XY is treated in a more general framework of the 2D Coulomb gas, 
having two kinds of quantum numbers for 'charges' and 'magnetic monopoles'|^ 
0. Several models, such as p-clock models, Ising model, three and four-state 
Potts models, Ashkin- Teller model are mapped in a unified way to the 2D 
Coulomb gas|ll[. 

Kadanoffp^ and Kadanoff and Brown[|l3j identified correlation functions of 
Gaussian, 8-vertex, and Ashkin- Teller models, whose critical dimensions vary 
continuously on the critical line. In the latter two models, only the few core- 
lation functions are known except at the decoupling point. They first made a 
comparison of the correlation functions of the three models at a special point 
on each critical line. Then, they used the marginal operator and the operator 
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product expansion to extend the stated connections on the whole critical line. 

The equivalence of the sine-Gordon model to the 2D Coulomb gas model has 
been shown by several authors Coleman showed the equivalence of the 

massive one-component Thirring model and the sine-Gordon model, order by 
order in perturbation expansion, and also showed the renormalizability. But his 
discussion failed in the region of 0^ > Stt. Luther and Emery[0 and Banks et 
al.|p^ showed the equivalence between the SU(2) massless Thirring model and 
the theory of bosons consisting of a free field plus a = Stt sine-Gordon model, 
which corresponds to the BKT line. So there is a hidden SU(2) symmetry at 
the BKT transition. 

A more systematic field theory treatment of the rcnormalization group cal- 
culation for the sine-Gordon model was performed by Amit et al.||l^. They 
resolved the problem met by Coleman, including the rcnormalization of the 
wave function. They calculated the higher terms beyond those of Kosterlitz 
and found the new universal quantity. 

The logarithmic corrections of the k = 1 SU(2) Wess-Zumino-Witten model, 
which is shown equivalent to the (3'^ — Stt sine-Gordon model, were system- 
atically studied by Affleck et al.^^. They foimd the universal relation of the 
ratios of logarithmic corrections to scaling amplitudes. This relation was used 
by Ziman and Schulz|21 for the problem of S=3/2 quantum Heisenberg chain. 

But this SU(2) symmetry is not apparent in the sine-Gordon model itself, 
and except for the BKT line, the symmetry is broken to 0(2) x Z2. How does 
the sine-Gordon model acquire an SU(2) symmetry on the BKT line? This 
problem, including logarithmic corrections, has been first treated by Giamarchi 
and Schulzp2|. They calculated the renormalized correlation functions, and 
found that the SU(2) symmetry of them is recovered on the BKT line. In their 
case the original model is apparently SU(2) symmetric on the (3'^ — Stt fixed 
point. 

There are several models mapped on the sine-Gordon model. Although the 
mappings are qualitatively correct, coupling constants and cut-off are renor- 
malized, so in order to determine the phase diagram one should use numerical 
calculations. But near the BKT transition, the divergence of the correlation 
length is essentially singular and also the logarithmic corrections exist. There- 
fore it is very difficult to find the critical point of the BKT type transition. In 
our previous paper using the level crossing of the eigenvalues of the transfer 
matrix or the corresponding quantum sine-Gordon Hamiltonian in ID, we have 
easily determined the transition point and identified the universality class. It 
was based on the SU(2) symmetry on the BKT transition line. 

In this paper, we perform the renormalization group calculation of corre- 
lation functions, whose critical dimensions become marginal on the BKT line. 
In the case that the SU(2) symmetry appears on BKT line, the 9 eigenvalues 
split as 5, 3 and 1-fold degenerate, i.e., the SU(2) multiplets structure. In other 
case, the 5 eigenvalues split as 3, 1 and 1-fold degenerate, and it is also possible 
to determine the BKT transition line by the level crossing of the eigenvalues. 
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In addition, new universal quantities are found. These may deepen our under- 
standing of the relation between the Abelian and the non-Abelian bosonization. 

This paper is organized as follows. In §2 the model is introduced, and the 
symmetry structure is discussed. In §3, we overview the correlation functions 
which become marginal at yo — y4, — 0. In §4, the renormalization equations 
are obtained for these functions. The hybridization between the marginal field 
and cos-\/8<?!> field is important. In §5, we consider the eigenvalue structure of 
the transfer matrix, and briefly summarize the results of §4. In §6, our results 
are applied to ID quantum and 2D classical systems. Section 7 is a conclusion. 



2 Sine-Gordon model 

The description of the symmetry and correlation functions of the Gaussian 



model in this and next sections is based on |23, E4, E5|. We first consider 



the 2D Gaussian model defined as the following Lagrangian 

= (2) 

The two-point correlation function for is 

2((/)(ri)(^(r2)) =-i^Relog(zi2/a), (3) 

where a is a short-distance (ultraviolet) cut-off and z = x + iy^ zi2 = zi — Z2. 
Strictly speaking, it should be also introduced a small mass /i to serve as an 
infrared cut-off, such as eq.(2.2) in |]l9j. The logarithmic behavior of the </> 
correlation function shows that it cannot be directly interpreted as a physical 
object. But the exponential operators of <j) behave as 

(exp(ie0(ri)) exp(-ze0(r2))> = |zi2/a|-^'^/2, (4) 

so they are candidates for the correlation functions of the critical theory. A 
convention needs to be explained regarding this formula. We did not include the 
divergent "self energy" factors coming from the terms in the exponent where 
the Green function is to be evaluated at 0. This means that we have really 
evaluated correlators of the "normal ordered" exponentials : exp(ie0) :. 

By the way, this Lagrangian is invariant under (/> — > -I- const, and 
—(j). This may be used to restrict configurations, with the identification of 
= 0-1- 27r/V2. In other words, one may compactify on a circle. In this case, 
the charges e are quantized as e = •\/2'T-(n:integer). One may also introduce 
new scaling fields exp(zTO-\/20(a;)), which create a discontinuity of by 27rm/-\/2 
around the point x. The two-point functions of them are 

2{e{r^)e{r2)) = --|Relog(zi2/a), 

2(0(ri)f?(r2)) = -^Imlog(zl2/a), (5) 
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and 



- -dy{iK9), dy4> = d,{iK0). 



(6) 



The field 9 is called a dual field to (p. This model is invariant under the transfor- 
mations 4> ~* (j)+const., 9 —^ 9 + const., that means U(l) x U(l) symmetry. The 
full symmetry group is extended to 0(2) x 0(2) by the discrete Z2 symmetries 
(z, (j), 9) (z, —(/), —9) and (z, 0, 9) (z, cj), —9). And there is a dual transfor- 
mation K ^ l/K, (p ^ 9, which exchanges the roles of electric and magnetic 
excitations. The self dual point = 1 is nothing but the k = 1 SU(2) x SU(2) 
WZW model. Another point of view is given by a next chiral decomposition 



1 



K 



K9 



1 



K 



'K9 



(7) 



This system is also chiral invariant. 

In this stage, it is a natural extension of the Gaussian model to introduce the 
interaction term cos ^/24>. Instead of this, in order to see the SU(2) symmetry 
explicitly on the BKT line, we consider the next sine-Gordon Lagrangian: 



C 



2ttK 



(V0)^ 



y4> 

27ra2 



(8) 



Note that the U(l) symmetry of cj) is explicitly broken to the discrete one </> — > 
Since this transformation divides the internal circle of (j) into two, 
the symmetry of this model is 0(2) x Z2 x Z2. The interaction term breaks the 
chiral symmetry. Furthermore this is invariant under — > (p + n/^/S, — s- —y^. 
However, calculations in this paper apply also to the usual BKT transition, 
simply by requiring the periodicity (f> = (j) + 2tt/^/8. In the latter case, the 
symmetry is 0(2) x Z2. 

From the perturbation expansion of the partition function, the renormaliza- 
tion group equations for the sine-Gordon model under a change of the cutoff 
a —> e'a up to the lowest order in yo and y^ are 



dl 

dy<i>{l) 
dl 



-yAi)yo{i), 



(9) 



where K — 1 + ij/o- For the finite system, I is related to L by e' = L. Note that 
there are three critical lines: y^ = {yo < 0) corresponding to the Gaussian 
fixed line, y^ = ztya {yo > 0) to the BKT lines. 
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3 Correlation functions on the Gaussian Fixed 
line 

We review the correlation functions on the Gaussian fixed \me{y^ — 0). The 
correlation functions of the Gaussian model are in general, 



{Ofi^rn (^1 —rL,~m 
,2 



exp 



- ( n^/C + ^ ) log(^) - 2mm(Arg(ri2) + J] 



On,m = exp(znV2(/>) exp(ir7iV20), 



(10) 



where Arg(ri2) = Imlog(zi2) is the angle of the ri — r2 vector. Thus On,m has 
a critical dimension Xn,m and a spin ln,m given by 



(11) 



The expectation value {TljOnj,mji'''j)) is zero unless the charge and monopole 
neutrality condition J2 "-j = 0; X] '^j = is satisfied llj, |l^, which reflects the 
underlying U(l) x U(l) symmetry. 

Then the flelds which become marginal (x = 2) and spin zero at t/o = 0{K = 
1) are the four fields of the form On.m , 

A={O2,0 + O-2.0) ^ \/2cos Vs./., ^(02,0 - O-2.0) = V2sin V8(/), 
v2 V 2 j 

O0.2, Oo,-2, (12) 

and the four descendant fields, which are expressed as the derivatives of the 
0±i ±1 fields. 



a[aOi,i +aO_i,i], a[50_i,_i +50i,_i]. 



(13) 



where d = dz, and we take the symmetrized and the antisymmctrizcd forms of 
the function of 6. There is one more field, that is 



M = 



K 



dx J \dy 



(14) 



The correlation functions for On,m are already shown. In order to calculate 
the correlation functions for descendant fields, it is enough to know 



^(90i4(ri)aO_i,_i(r2)) 



exp 



vl 



log(ri2/a) 



(15) 
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where — A{K + — 2) ~ ?/q. To calculate the correlation function for the M. 
field, we introduce the next fields 



a 
K 



dx dy 



a 



2a(9)</), 



(16) 



then 



(7^+(rl)7^+(r2)> 
(7^+(rl)7^-(r2)> 



fl2 

a 



, (7^-(rl)7^-(r2)> 



fl2 

a 



0. 



(17) 



That is, are the fields with the critical dimension 1 and the spin ±1. There- 
fore, 

{Min)M{r2)) - (7^+(rl)7^-(rl)7^+(r2)7^-(r2)) = cxp[-41og(ri2/a)]. (18) 

that is, the critical dimension of this field is always marginal(x — 2). Here we 
have implicitly taken the normal ordered form Ai =:TZ'^TZ~ : for the marginal 
field. The M is nothing but the marginal field of Kadanoff et al.[^ |l^, and 
7^+*^^^ correspond to their -F'i,o(-fo,i)- Note that in these nine fields, although the 
critical dimensions of On,m and their descendant fields vary with the parameter 
2/0, the marginal field M has always the critical dimension x = 2. The next 
equations 

(7^>r)^(.2)) = -^(^)"^(7^->^o 

are useful for further calculations. The three point functions are 



(19) 



{M{n)M{r2)Mir3)) = 0, {M{n)M{r2)On,M) = 



(20) 



The former comes from the Wick's theorem and (|17|), (|19|), and the latter from 
the neutrality condition. 

In the usual BKT transition where the SU(2) symmetry is not explicit, the 
difference is only that the critical dimension of the descendant fields is 2 + 1 + 
1/8 = 25/8 at 2/0 = Uip — 0- Therefore only the five fields are marginal at 
the multi-critical point. The usual BKT transition point can be obtained by 
modding out the SU(2) symmetric case with a Z2 symmetry[23, |4, 25 . 



4 Renormalization of correlation functions 

In this section, we proceed to include the interaction y^ cos \/80. When we 
calculate the correlation functions, there appear divergences coming from both 
a short range and a long range, because of the nature of the 2D Green function. 
To treat such a problem, we use the renormalization treatment. The correlation 
functions for (exp(i-v/86'i) exp(— «\/802)) and (\/2 sin(\/8(/)i)\/2 sin(V8(/)2)) (for 
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brevity we will use (jji for 4>{ri)) have been already obtained by Giamarchi and 
Schulz|2^. Their results are 



i?2 



2(sin(\/8(/)i) sin(\/8(^2))/ = C2 exp 



ln(r/a) 



i?3 = (exp(i'\/80i) exp(— iV8(?2))/ = C3 exp 
i?4 = (exp(-i\/86'i)exp(iV86'2)>/ = i?3. 



dl{A + 2yoil)) 

\n{r / a) 

dl{4~2yo{l)) 



(21) 



where Cj are the integration constants depending on the regularization method. 

4.1 Correlation functions of the marginal and cos ^/Scp fields 

As for the fields Ai and \/2 cos y/Scf), they are hybridized by the interaction term 
as 



_M2 cos 

= f ^(A^2(<^i-</'3)')(exp(zV80i)expHV803)) 

= -Kj -^oc\z^^ - ^3"2')(%^2' - ^3"2')exp(-4i^log(r3i/a)). (22) 

To derive this results, we use Wick's theorem [pO[, the U(l) symmetry and (|l9|). 
Terms which are not invariant under the global transformation + const. 

should be zero. The divergent parts relating to renormalization appear near 
^si ri2 and r32 < ri2. When r^i <^ ri2, z^^ — ^i^ii^ ~ ^3i/-Zi2), therefore 
the integrand of (E2|) is 



Zl2 




Z31 


a 




a 



exp{-4K log(r3i/a)). 



When r32 ri2, the divergent part of the integrand is 



(23) 



-232 



exp(-4i^ log(ri2/Q;)) 



^12 


-4 


^32 


a 




a 



exp(-2yolog(ri2/a)), (24) 



where we use K — 1 + ya/2. Note that only the terms which contain |2;32|~^ 
contribute to the divergence of the integral; other terms such as z^2 cancel 
by the integral. In order to treat these divergences, we exclude two circles of 
radius a around ri and r2 from the domain of integration over . Then, with 
the change of cutoff a' — ae'" , this integral is renormalized as 



^exp(41og(ri2/a)) / 



_M2 cos 
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= ^ exp(41og(ri2/a')) / ^^-^(cos%/8(?!)iA^2 cos^ 
27r Ja' a'-^ 

- 2yl[l-yo\og{n2/a)]dl. (25) 
4.1.1 Near the Gaussian fixed line 

At first we treat the case of Ij/^/yol ^ 1- Let us consider tlie two liybridized 
states between the marginal and cos \/8(/) fields, 

A = M+ a{y^/yo)V2cos VSc/), 

B = V2cosV8(l) + b{y^/yo)M. (26) 

the orthogonal condition We consider the condition that the correlation 
function {A1B2) stays zero under renormalization. 

{AiB2)i = (y0/2/o)exp(-41og(ri2/a))[a(l - 2?/olog(ri2/Q;)) + fe] 

f ^V2[(cos V8^iM2 cos + abiy^/yof{l ^ 2)] 

+higher order terms. (27) 

By using y^{l) ~ y^{0) exp{—yo{0)l) and yo{l) — 2/o(0), the function defined as 
F = (AiB2)/exp(41og(ri2/a))exp(yo(O)/)(2/o(O)/2/0(O)) behaves approximately 
constant a + b for small enough yo,y^. In order to set {A1B2) = 0, first of all 
it should be a + 6 = 0. In addition, in the course of the renormalization, other 
terms may appear. For the infinitesimal transformation a' = ae'^\ using eq. 
we obtain 

F = F' + [-202/0 + (2/o/y0)2\/2y^(l + ab{y^/ (j)of)]dl, (28) 

where F' is the function F with the new value of a. Remark that in the course of 
the renormalization of j/o ^C)g{ri2/ a), the term y^ log(ri2/Q!)(iZ appears, however 
this disappears by the cancellation with the term in higher order expansions 
^ . Then the necessary conditions for F = Q under renormalization are 

a + b = 0, 

-2ayo + 2V22/o(l + a^'(y0/0o)') = 0. (29) 
The solution for them is 

a = -b = ^+0{{y^/yof). (30) 
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renormalized correlation function for the marginal-like field Correla- 
tion function of the marginal A^-like field is 

= exp(~41og(ri2/a))[l + a2(y^/y„)2(i„2y„log(ri2/a))] (31) 

/_j2 
^[(cos V8</)iM2 cos V8^3) + (1 ^ 2)]. 

The function defined as F = (^1^2)/ exp(41og(ri2/Q;)) behaves approximately 
constant 1. For the infinitesimal transformation a' = ae''', using eq. (psf), we 
obtain 

F = ex.^{[~2a?{y^/yfifyo + 2a?{y^/yofyQ + 4V2a(?/0/yo)y0)]dO x F' 
= exp{8yo{y4,/yofdl)xF'. (32) 

As a result, 

^ln(r/ q) 



Ro = (A1A2)/ = Coexp 



dl[A-8yo{l)iy^/yQ) 



(33) 



renormalized correlation function for the cos v 8(/i-like field Correlation 
function of the cos \/8(h-\ikc field is 



(BiB2)7 = exp(-41og(ri2/a))[l-2yolog(ri2/a) + 62(y^/2/o)'] (34) 

, , ., ,(cosV80i>l2Cos%/8(?!)3) + (1 ^ 2)]. 
2tt J 

The function defined as = (i?ii?2)/ exp(41og(ri2/a)) behaves approximately 
constant 1. For the infinitesimal transformation a' = ae'*', 



F = cxp([-2t/o + (?y0/2;o)'yo(252 + 4A/26)]dO xF' 
= exp([-2?/o - 4.yo{y4,/yof]dl) x F'. 



As a result, 

i?i = (B1B2)/ = Ciexp 



ln(r/ q) 



dZ[4 + 2yo(O(l + 2(2/0/yo)')] 



(35) 



(36) 



4.1.2 Near the BKT transition line 

Next we treat near the BKT transition, that is, = ±?/o(l + t), \t\ <C 1. With 
this parametrization, t plays the role of the deviation from the critical point, 
such as (T — Tc)/Tc . Let us consider the two combination of the marginal and 
cos VS^ fields, 

A = M + a\/2 cos \/8(f>, 

B = V2cosV8(j) + bM. (37) 
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the orthogonal condition The correlation function {AiB2)i is obtained by 
replacing y,j,/yo = 1 in the previous subsection. The function defined as F = 

(Aii?2)/ cxp(41og(ri2/Q!)) behaves approximately constant. The conditions for 
F — under renormalization are 



a + b = 0, 
^2ayo + 2V2y4,{l + ab) = 0. 



(38) 



The solution for them is 



a = -b= ±— for y^ = ±yo{l + 1). 



(39) 



renormalized correlation function for the marginal-like field The func- 
tion defined as F = (A1A2)/ exp(41og(ri2/Q:))/(l + b"^) behaves approximately 
constant 1. For the infinitesimal transformation a' = ae'^K 



exp 



= exp 



-2a'^yo + 4V2a?/0 



1 + 62 



dl\ xF' 



2yo(l + -t] dl 



X F'. 



As a result, 

i?o = {AiA2)i = Co exp 



-L 



\n.{r / a) 



dl 



4-22/0(01 1+3* 



(40) 



(41) 



renormalized correlation function for the cos \/8(/)-like field The func- 
tion defined as F = (i?ii32)/ exp(41og(ri2/a))/(l -|- b^) behaves approximately 
constant 1. For the infinitesimal transformation a' = ae'^\ 



F = 



exp 



-2yo + 4V26y0 



= exp 

As a results, 

i?i = (BiB2)/ = Coexp 



1 + 62 

-4yo [l+'^t\dl 



dl] xF' 



X F.' 



ln(r/a) 



dl 



4 + 4t/o(0( 1+3* 



(42) 



(43) 
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4.2 Correlation functions of the descendant fields 

The renormalization calculations for descendant fields {dOi.i etc.) are straight- 
forward extension of the method by Giamarchi and Schulz[2^ 

a^idOis + 90-i4)(^i)(50i,_i + aO_i,_i)(r2))/ 

(fx3 



y± 

27r 



^(90i,i(ri)aOi,-i(r2)02,o(r3)). 



(44) 



The first two terms are already calculated. Therefore it is enough to estimate 
the divergent part of the integrand of the third term, 



a2(aOi,i(ri)50i,_i(r2)02,o(r3)) 



^31 

a 



X exp 



-2 ^Z32^-2 



1 yQ\ 1 a 

2 4 zi2 2 Z31 



2 4 zi2 



1 a 

2 2:32 



T^i + ivo- ^) log(ri2/a) - yolog(r3i/a) - yo log(r32/Q;) 



(45) 



This diverges near r^i ^ ri2 and r32 ^ ri2. In the case of r^i ^ ri2, 
Zi2i^ — Z31/Z12), thus the divergent part is 



Vo 
4 



Z31 


-2 


^12 


-4 








exp 


a 




a 





-^log(ri2/a) 



(46) 



Note that only the terms which contain jzisj"^ contribute to the divergence. 
The divergent part near r23 <C ri2 is the same form. 
Let us consider the next function 



F = ((aoi,i + ao_i,i)(ri)(aoi,_i + ao_i^_i)(r2))/ 

-1 

exp 



X ( 2^ + ^ 



)log(ri2/a) 



(47) 



which reduces the constant 1 when y,p — yo ~ Q- Then considering renormaliza- 
tion behavior, we obtain as before 



As a results, 

Rr 



F = exp{-2y^dl) x F' . 

((aoi4 + ao-i,i)(ri)(aoi,-i + ao_i,_i)(r2))/ 

/•log(ri2/a) 

Csexp - / {4 + 2y4l))dl , 

((aoi,_i + ao_i,_i)(ri)(aoi.i + ao_i.i)(r2))/ 



(48) 



(49) 



i?5 
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The calculation is the same for the antisymmetric combination, except for the 
sign of 2/0, 



i?7 = 

= C-j exp 
^8 — Rr- 



((aOi,i - 90_i,i)(ri)(50i,_i - aO_i,_i)(r2))/ 

log(r-i2/Q) 

{A-2y^{l))dl , 



(50) 



5 Eigenvalue structure 

Conformal field theory ^ is an efficient method to determine critical di- 
mensions of 2D systems. One of the most useful application of this theory is for 
the finite-size scafing. When we denote the transfer matrix of a strip of width 
L, with periodic boundary condition, by exp(— iJ), then the eigenvalues En are 
related to the scaling dimension a;„ of the scafing operators of the theory as[p9t 

2'HXn 

in the fimit of L — )■ cx). But this relation is exact only at the fixed point. In 
general, there appear corrections by terms involving irrelevant (marginal) fields. 

It is possible to relate the eigenvalues of the transfer matrix with the renor- 
malized correlation functions obtained in the previous section. The renormalized 
critical exponents 77,1 (/) are related to the correlation functions aspQ] 



E^{L) - E,{L) = ^ (51) 



Rn = exp 



ln(r/a) 

dlr]n{l) 







(52) 



and by the use of the eq.(|51|), which relates the critical exponents 77 to the 
eigenvalues for the finite size system, we obtain 

= Xn{l) = ^Vn{l)- (53) 

Although the eq.(|5l|) is satisfied under the condition of scale invariance, by the 
use of the renormalization group, we extend this relation to the region where 
scale invariance is not strictly satisfied but the system size is sufficiently smaller 
than the correlation length ^ of the L — > 00 limit. The renormalized critical 
dimensions are (close to the BKT transition) 

Xo{l) - 2-yo{l){l + ^t), 

xiil) - 2 + 2yo{l)il + h), 
X2{1) = 2 + yo(0, 

x^{l) = x^il) = 2-yoil), (54) 
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and 



xr{l) = xsil)^2-y^{l)^2TyQil + t). (55) 

These results mean, on the BKT transition Une (for example ytti{l) = 2/o(0) 
that the eigenvalues of the transfer matrix corresponding to the fields xo{l), X3{1), Xi 
and X'!{l),xs{l) become degenerate, also those to 0:2(0 ^'^^ xz{1),xq{1) degen- 
erate. This structure reflects the fact that the 0^ = Stt sine-Gordon model 
corresponds to the SU(2) massless Thirring model||l8| or the SU(2) k — \ Wess- 
Zumino-Witten modellg^l- On the BKT line, yo(0 renormalized as yo(0 — 
1/logL. The ratios of the logarithmic correction terms in x^{l),X2{l),Xi{l) are 
-1 : 1 : 2, in agreement with the SU(2) k = l WZW model[|o|. Although the 
convergence of the logarithmic term is very slow, this relation can be used to 
eliminate the logarithmic correction |^ p6| . 

In the neighborhood of the BKT transition line, in the correction terms 
of xq{1),x\{1) and xz{1),xq{1),xt{1),xq{1) there appear the terms linear to the 
distance t from the BKT line, and the ratios of them are — | : | : 1 : — 1, 
indicating the new universal relations. Moreover, for example, xq{1) — x^{l) is 
linear to t, a useful relation to determine the BKT critical line. This describes 
how the SU(2) symmetry breaks down to the U(l) x Z2 symmetry in the Abclian 
bosonization. 

Close to the Gaussian fixed line, the only differences are 
xo{l) = 2-4yo(O(y0(O/yo(O)', 

xi(0 = 2 + yo(O[l + 2(j/0(O/yo(O)'], (56) 

therefore the ratios of the coefficients (ycti/yo)^ terms are —4 : 2. In this case, the 
difference 0:5 (Z) — xt{1) is linear of the deviation 2/^(0 from the Gaussian fixed 
line, since they exchange each other under the transformation — > (/) -I- tt/^/S. 
And the difference xi{l) — 2:2 (0 is quadratic of y^. This relation may be used 
to determine the Gaussian fixed line. 

Finally, a comment on the symmetry. The SU(2) symmetry is broken to 
0(2) X Z2, while the symmetry structure of the model (||) is 0(2) x Z2 x Z2, 
therefore an additional Z2 symmetry. This additional Z2 symmetry in SU(2) x Z2 
on the BKT line is needed to inhibit the SU(2) symmetric relevant fieldlSg]. In 
general, for the SU(n) critical model, an additional Zn symmetry is necessary 
to stabilize the massless phase . 

In the usual BKT transition, where SU(2) symmetry is not explicit, the 
difference is only that the critical dimension of the descendant fields becomes 
25/8, so they are not marginal. Nevertheless, there remains degeneracy includ- 
ing logarithmic term between the exp{±i^/89) fields and the marginal-like field 
on the BKT line. This fact is useful to determine the BKT transition point 
from eigenvalues. 
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6 Physical systems 



In this section, we treat an S=l/2 XXZ chain with next-nearest neighbor interac- 
tions which corresponds to the sine-Gordon model with an 0(2) x Z2 x Z2 sym- 
metry. As an example of the sine-Gordon model with a simple 0(2) x Z2 sym- 
metry, we treat a bond-alternation S=l/2 chain. And we treat the 2D classical 
p-clock model as an example of the sine-Gordon model with an 0(2) x Z2 x Zp 
symmetry. 

6.1 S=l/2 XXZ chain with next-nearest-neighbor inter- 
action 

In the previous workp^, we have studied the S = 1/2 XXZ chain with a 
competing-interaction: 

L L 

H = $^(5J5^i + Sj^jVi + A5|5|+i) + a^(5J5j+2 + S]Sy^, + AS^S^^,), 

(57) 

with the periodic boundary condition Sl+i = Si, and L is the number of spins 
{L even). 

The above model is exactly solved on lines a = and a = i. On the line 
— 1<A<1,q; = 0, the ground state is the spin-fluid state, characterized 
with the gapless excitation and power-law decay of correlation functions. In 
the region A > 1, a = 0, this system is Neel ordered and it has a two- fold 



degenerate ground state with an energy gap|31, |3^. On the line a = ^, the 
ground state is purely dimerized [^3[ It is also proven the existence of the 
energy gap and the uniqueness of the two- fold degenerate ground state The 
dimer state is characterized by the excitation gap, the exponential decay of the 
spin correlation, and the dimer long-range order. 

Let us examine the symmetry of the Hamiltonian(|57|) . This model is in- 
variant under spin rotation around the z-axis, translation(5j'^'") ^ S'^iY^), 

space inversion(S'J^^'^^ S^L-f+i)^ spin reverse(5'| —Sj.Sf ~^f)^ 
and conjugate(S'| — + Sj,S^ S^). Therefore eigenstates are character- 
ized by z-component of the total spin(S'|i = X^'^'J): wavenumber(q = 2t:/L), 
parity(P — ±1), spin reverse(T = ±1), and charge conjugate C. The charge 
conjugate C is redundant, because of CPT = 1 as will be shown later. For 
L = 4n, the ground state is a singlet(S'|, = 0,q = 0, P = 1,T = 1). The symme- 
tries of several low excitations are classified in Table 1. The operators in spin 
representation are also shown. 

Next we consider a corresponding sine- Gordon model. After a Jordan- 
Wigner transformation, the model(p7|) is transformed to the ID spinless fermion 
system. The continuum limit of this is a Tomonaga-Luttinger liquid p^, or 
equivalently a sine-Gordon modelp^. Using the same procedure, we obtain 



15 



expressions in the sine-Gordon model which correspond to the spin opera- 
tors(Table 1). The marginal-like field and the cos\/80-like field are parts of 
the Lagrangian, so they have the same symmetry with the ground state, and 
the corresponding spin operators are the parts of the Hamiltonian with the 
same symmetry. Except for the Gaussian fixed line, a hybridization occurs 
between the marginal and the cos \/80 fields. Note that the fields in the sine- 
Gordon model are defined on the infinite plane, whereas operators in the spin 
model are defined on the cylinder. The former can be mapped to the latter 
by fiz) = L/27rlog(z). The symmetry operation in the sine-Gordon model 
corresponding to the spin reverse (T) is 

''(j) + TT/V2, 9 ^ -e + TT/V2, (58) 

the operation to the space inversion (P) is 

(1)^ -(l) + Tr/V2, 9 ^ e + Tr/V2, z ^ z, (59) 

and the operation to the charge conjugate (C) is 

(j) ^ 4>, 9 ^0, z ^ z, (60) 

therefore the successive transformation of them is the identity CPT =1. In 
addition, the operation corresponding to the translation by one site is 

(j)^(l3 + TT/V2, 9^9 + 7r/V2. (61) 

The symmetry breaking related to the phase transition is as follows. The 
sine-Gordon model is invariant under (f) — > 0-1- 27r/'\/8, which means expectation 
values of {ainy/2(f>), (cos \/20) are zero in a symmetrical phase. In the spin- fluid 
region, the is renormalized to zero, so no symmetry breaking occurs. In the 
dimer region, — > +oo, therefore one has a long-rang order of the field, 
whereas correlations of 6 decay exponentially. The average value of the ordered 
field is (0) = tt/\/8. In the spin system, this corresponds to the symmetry 
breaking of the translation invariance. In the Neel region, y^ — > —oo, ((/>) = 0. 
This corresponds to the symmetry breaking of the translation invariance, space 
inversion and spin reverse. 

We compare the renormalization calculations in sections 4,5 with the nu- 
merical results. The whole phase diagram is shown in Fig.l. The normalized 
excitations LAE/2ttv for A = 0.5 are shown in Fig. 2. The general behavior 
of them are consistent with the renormalization calculation. It is seen that the 
SU(2) symmetry appears on the BKT line|26|. Secondly we investigate the ra- 
tios of logarithmic terms. From eqs. (|5^), on the BKT line {t — 0), by taking 
averages 

^[x,{l) + x,{l)], i[xi(0 + 2x3(0], (62) 
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we can eliminate the contribution of the logarithmic corrections, and simulta- 
neously we can confirm the ratios of them. In Fig. 3, these averages are shown 
as a function of a system size L. As expected, they converge to 2 with 
due to the irrelevant field L-2L-2I {x = 4)|^. The extrapolated values are 2 
within 0.2% errors, comparing the bare values of x„(Z)(5-15% errors). Finally, 
we examine the ratios of the linear terms of t in the critical dimensions. From 
eqs. (jsj), (|55|), by taking averages 

^[xoil) + xi{l) + xsil)], i[x5(/)+X7(Z)], l[lxo{l) + X5{l) + \x3{l)], (63) 

the linear terms of t should be annihilated. In fact, as is shown in Fig. 4, in the 
neighborhood of the critical point ac = 0.2764, the linear components of t are 
almost absent(the points shown in Fig. 3 have been extrapolated as l/L"^). The 
coefficients of t linear terms are at least 10~^ less than those in the raw data 

Xn{l)- 

6.2 Bond-alternation S=l/2 spin chain 

This model is described by the following Hamiltonian 

^ = E(l + 5{-iy){Sp^^, + Sj^jVi + A5|5|+i). (64) 

After a bosonization, we obtain 

The correspondence between the spin variable and the sine-Gordon model is 
the same as before (Table 1). The symmetry structure in this model is a simple 
0(2) X Z2. At A = —l/y/2,5 = 0, the critical dimension of sin\/2</> becomes 
marginal. Therefore, in the neighborhood of this point, there appears the BKT 
transition[Q. Although there exist higher terms such as cos\/80, we neglect 
them for simplicity. 

After a simple transformation of the fields (j) 2(f>, 9 9/2, we can use 
the results of the sections 4 and 5, except that the critical dimension of the 
descendant fields becomes x ~ 25/8 at yo = = 0. The corresponding spin 
operators which become marginal at yo — y^, — are {~iy Sj, (— 1)-' (5^5*^1^^ + 
S~S+^^), S+ S+^j^S+^^S+^^iS^ = 4,g = 0,P = 1) and the term corresponding 
to the marginal field. 

In the spin-fluid region, the y^ is renormalized to zero, so no symmetry 
breaking occurs. As for the region where y^ flows to inflnity, one has long-range 
order of the field (0) = =F7r/V8 for y^ — > ±00. This long range order is 
not related with any symmetry breaking. By using the symmetry consideration 
of the previous subsection, this corresponds to the symmetry breaking of the 
translation by one site. However, the translation invariance is already broken 
in the original model itself, there is no symmetry breaking in the whole region. 
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6.3 p-clock model 

We consider the case of the 0(2) x Z2 x Zp symmetry sine-Gordon model 

The symmetry of (p is invariant under — > + 2it/p^/2, i.e., sifting the circle 
coordinate </> by \/p times its period(27r/-\/2). The cosp-\/2(^ term becomes 
marginal at = 0, if = By parametrizing Kp^/A = 1 + yo/2, the 

renormahzation equations are the same as (^. The marginal fields at j/o = 
2/0 = are the three, cospv^^, sm.p\/2(j), A4. The renormahzation behavior 
of them has been described in §4,5. There is no degeneracy on the BKT line. 
However, the renormalized critical dimension for Oo.m | |2^ is 

2:^(0 = ^(1-^2/0(0), (67) 

therefore the ratios of the xo{l) and x'j^{l) is rational 2 : including 
logarithmic corrections. It is possible to use this to determine the BKT critical 
line, because xo{l) — lQx'^{l)/m?p^ is linear to the deviation t from the BKT 
line. The ratios of the logarithmic corrections of them may support to determine 
the critical dimensions and to check the consistency. 

When p is even, besides the above one, another relation appears. The renor- 
malized critical dimension for cos p4>/\/2 isp2[ 

i(l + |2;o(0(l + |0) fort/0-yo(l + t) 
\{l-\vo(l){l + 2t)) fory0--yo(l + i), ^ ^ 

and the renormalized critical dimension for sinp(j)/^/2 is 

^'m-J |(l"|2/o(0(l + 2t)) for2/0 = yo(l + t) ..q^ 

Therefore the lower part of them can be used to determine the BKT line as 
xo(0-4x'(0. 

What are the corresponding real systems? One candidate is the p-clock 
modcl[|l^, which is equivalent to the 2D XY model with Zp anisotropy. In this 
model the spins at each site can take only p-discrete angles 27rZ/p, I — 1, . . . ,p. 
The classical Hamiltonian is 



H = —K cos 

{r,r') 



27r 



^ p 



m-nP) , (70) 
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where the sums over r index the sites of a 2D lattice, the symbol (r, r') indicates 
a sum over nearest neighbor lattice sites only, and K = J/ksT. Note that this 
is a classical Hamiltonian, so it should be interpreted as a Lagrangian, rather 
than the logarithm of the transfer matrix. 

This model has two transition points. The upper critical point is between 
the disordered and the 2D XY like phase, and the lower critical point is from the 
2D XY like phase to the Zp symmetry breaking phase. The upper critical point 
is described by the usual O2 x Z2 BKT transition, whereas the universality class 
of the lower critical point is the O2 x Z2 x Zp BKT transition described in this 
subsection. 



7 Conclusions: level spectroscopy 

The idea that the level crossings of the low excitations can be used to deter- 
mine the critical point (hereafter we call it as "level spectroscopy") originates 
from the work by Ziman and Schulz[pT|, who studied the S=3/2 isotropic spin 
chain on the basis of the conformal field theory and the renormalization group. 
The problem how the system becomes SU(2) symmetric on BKT line from the 
anisotropic phase was treated by Giamarchi and Schulzp^, and this was used 
for the problem of the S=l/2 NNNI chain[^. We have proceeded in this pa- 
per in the case that the SU(2) symmetry is not apparent on the BKT line, by 
considering the hybridization between the marginal and cos fields. 

The level spectroscopy method is completely different from the finite-size 
scaling |4^, and for the BKT problem, the former is superior than the latter. 
In the finite-size scaling, one uses the data from several lengths to construct 
a scaling flow relation, then one searches the fixed point. But in the BKT 
problem, there are continuous fixed points below the critical temperature, so 
it is difficult to determine the BKT critical point by the finite-size scaling. In 
addition, there is the problem that the correlation length diverges singularly, 
and logarithmic corrections. 

On the contrary, in the level spectroscopy, it is used the symmetric struc- 
ture of the eigenvalues of the transfer matrix in order to determine the critical 
point and to obtain the critical dimensions. The renormalization process is 
already done by hand, not by the numerical data. Also the singular behav- 
ior of correction terms can be eliminated. In this case, only the data in one 
length are needed. In general, there exist the corrections from the irrelevant 
field L-2L-2^{x = 4), therefore it is needed extrapolations. Nevertheless the 
convergences are extremely fast. 

Finally, although it is possible to consider the sine-Gordon model with the 
0(2) X Z2 X Zp symmetry, the SU(2) symmetry appears on the BKT line only 
the 0(2) X Z2 X Z2 case. Such a symmetry does not occur in other cases. The 
usual 0(2) X Z2 BKT critical point is also special, because in this point the 
Gaussian and the orbifold model are equivalent p3, M. 
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Table 1: Identification of eigenstate of XXZ model in the sine-Gordon language. 



Symmetries of Identification in Identification in 

eigenstate spin language sine-Gordon model 



<1 


iS'j-' 


T 


P 












1 


1 


1 


1 


TT 


1 


* 


-1 






TT 





-1 


-1 




Oi.o + 0-1,0 = 2 cos V2(j) 


TT 





1 


1 




i[Oi,o-0-i,o] =2sinV2</. 


2tt/L 





-1 


* 




del) 


2tt/L 


1 


* 


* 


exp(27rij/L)S't 


Oi,i 





2 


* 


1 




Oo,2 





1 


* 


1 












1 


1 


a part of the Hamiltonian 


M 





1 


* 


-1 




i[aoi,i-ao_i,i] 








-1 


-1 


— {Sj'Sj^i + Sj Sj[^i)Sj_^_2 


i[O2,0-O_2,0] 








1 


1 


a part of the Hamiltonian 


02,0 + 0-2fl 
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